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^) METHODE POUR DEFORIMER GRADUELLEMENT UN MODELE STOCHASTIQUE D'UN MILIEU 
HETEROGENE TEL QU'UNE ZONE SOUTERRAINE. 

- Methode pour deformer de fagon graduelle, en totar preservant la continuity entre ces parties, 
lite ou localement un modele stochastlque de type gaussien - Application a la construction de modeles stochastiques 
ou apparente d'un milieu heterogene tel qu'une zone sou- de reservoirs contraints ^ des donnees non linialres (don- 
ten-aine, contraint par un ensemble de parametres relatifs h nees liees a I'ecoulement des fluides). 
la structure du milieu. 

- La methode comporte le tirage d'un nombre p au moins 
egal a deux, de realisations independantes d'une partie au 
molns du modele choisi du milieu, parmi toutes les realisa- 
tions possibles et une combinaison lineaire de ces p realisa- 
tions avec p coefficients tels que la somme de leurs carr6s 
est 6gale k 1 . Cette combinaison lineaire constitue une nou- 
veile realisation du module stochastlque et elle se d^forme 
graduellement lorsque Ton modifie graduellement les p 
coefficients. La methode peut comporter plus gen^ralement 
plusieurs etapes iteratives de deformation graduelle, avec 
combinaison a chaque etape d'une realisation composite 
obtenue a Tetape precedente avec q nouvelles realisations 
independantes tiroes de I'ensemble des realisations. La me- 
thode donne la possibilite de deformer graduellement des 
realisations d'un modeie representatif du milieu tout en mo- 
difiant.les parametres statistiques relatifs a la structure du 
milieu. La methode permet aussi des deformations graduel- 
les individueiies de differentes parties du modeie tout en 
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La presente invention conceme une methode pour deformer 
graduellement partiellennent ou en totalite des realisations ou 
representations d'un milieu heterogene tel qu'une zone souterraine, qui sont 
toutes la traduction d'un modele stochastique de type gaussien ou apparente 
5 c'est-a-dire incluant un modele gaussien sous jacent. 

La methode selon Tinvention trouve des applications notamment dans 
la construction d'un modele stochastique d'une formation souterraine, 
contraint par des donnees non lineaires. 

ART ANTERIEUR 

10 Des exemples d'utilisation de modeles gaussiens pour modeliser la 

structure du sous-sol sont decrits par : 

- Joumel, A.G. and Huijbregts ChJ.: "Mining Geostatistics", Academic 
Press. 1 978 ; ou 

- Matheron, G. et a! : « Conditional Simulation of the Geometry of Fluvio- 
15 deltaic Resen/oirs », paper SPE 16753; SPE Annual Technical 

Conference and Exhibition, Las Vegas, 1987. 

La contrainte des modeles stochastiques a un ensemble de donnees 
non lineaires peut etre consideree comme un probleme d'optimisation avec 
definition d'une fonction objectif mesurant Taccord entre des donnees 

20 mesurees dans le milieu a modeliser et des reponses correspondantes du 
modele stochastique, et minimisation de cette fonction. Une bonne methode 
de deformation se doit de respecter les trois caracteristiques suivantes : a) le 
modele deforme doit etre une realisation du modele stochastique ; b) il doit 
produire une fonction objectif reguliere pouvant etre traitee par un algorithme 

25 d'optimisation efficace ; et c) la deformation doit etre telle que Ton puisse 
couvrir en entier Tespace des solutions. 
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Parmi les methodes de deformation connues, on peut citer la methode 
selon laquelle on echange les valeurs en deux points d'une realisation du 
modele stochastique i.e. d'une realisation ou representation du milieu etudie 
traduisant le modele. Cette methode est utilisee dans la technique dite de 
5 « recuit simule » bien connue des specialistes. Un exemple en est decrit 
par : 

- Deutch, C: "Conditioning resen/oir models to well test information", in 
Soares» A. (ed.), Geostatistics Troia'92, pp.505-518, Kluwer Academic 
Pub., 1993. 

10 Comme un echange arbitraire viole la continuity spatiale d'un modele 

stochastique. il est necessaire d'inclure un variogramme dans la fonction 
objectif, ce qui rend Toptimisation tres fastidieuse. 

Une autre methode connue est la methode dite des points pilotes, qui 
est decrite par exemple par : 

15 -de Marsily, G. et al : « Interpretation of Interference Tests in a Well Field 
using Qeostatistical Techniques to fit the Permeability Distribution in a 
Reservoir Model » in Verly, G. et al (ed.). Geostatistics for Natural 
Resources Charaterization, Part 2. 831-849, D. Reidel Pub. Co. 1984. 

Cette methode consiste essentiellement a choisir un certain nombre 
20 de points d'une realisation (representation) initiale. a calculer les derivees 
^jes coefficients de sensibilite) de la fonction objectif par rapport aux valeurs 
en ces points, et ensuite a modifier les valeurs aux points pilotes pour tenir 
compte de ces coefficients de sensibilite. Une nouvelle realisation est formee 
par la methode connue dite de krigeage conditionnel. II peut en resulter 
2^, qu'une modification des valeurs aux points pilotes en utilisant les coefficients 
de sensibilite. conduise a une modification indue du variogramme meme si le 
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krigeage conditionnel est utilise ensuite. en particuller quand le nombre de 
points pilotes devient important. 

DEFINITION DE UINVENTION 

La nnethode selon I'invention permet d'effectuer une deformation 
5 graduelle au moins partielle de realisation (ou representation) d'un modele 
stochastique de type gaussien ou apparente, d'un milieu heterogene tel 
qu'une zone souterraine, contraint par un ensemble de parametres relatifs a 
la structure du milieu. Elle est caracteris6e en ce qu'elle comporte le tirage 
d'un nombre p au moins egal a deux, de realisations independantes d'une 

10 partie au moins du modele choisi du milieu, parmi toutes les realisations 
possibles et une combinaison lineaire de ces p realisations avec p 
coefficients tels que la somme de leurs carres est egale a 1. Cette 
combinaison lineaire constitue une nouvelie realisation du modele 
stochastique et elle se deforme graduellement lorsque Ton modifie 

15 graduellement les p coefficients. 

L!a methode peut comporter plus generalement plusieurs etapes 
iteratives de deformation graduelle, avec combinaison a chaque etape d'une 
realisation composite obtenue a I'etape precedents avec q nouvelles 
realisations independantes tirees de Tensemble des realisations. 

20 La methode permet la deformation graduelle d'un modele 

representatif du milieu tout en modifiant les parametres statistiques relatifs a 
la structure du milieu. 

La methode permet aussi d'effectuer des deformations graduelles 
individuelles de differentes parties du modele tout en preservant la continuite 
25 entre ces parties (c'est-a-dire par exemple la reproduction d'un variogramme 
attache a ce modele). 
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De preference, les coefficients de la connbinaison lineaire sont des 
fonctions trigonometriques. 

Comme modele, on peut utiliser un modele gaussien ou apparente au 
type gaussien : modele lognormal. gaussien tronque, etc. 

Par rapport a Tart anterieur, la methode selon I'invention est fondee 
sur une base theorique plus solide pour la deformation des modeles 
stochastiques. Elle donne aussi la possibilite de modifier les parametres 
(souvent de valeurs incertaines a priori) de variogrammes en meme temps 
que la deformation d'une realisation d'un modele stochastique. Enfin, grace 
a la possibilite des deformations individuelles de differentes parties du 
modele, la methode selon I'invention donne une plus grande souplesse et 
une plus grande efficacite a I'operateur pour Tajustement du modele 
stochastique a la realite des terrains. En cela, la methode selon I'invention 
permet de rapprocher Tajustement des modeles stochastiques de 
Tajustement des modeles deterministes par zonation traditionnellement 
employe par les ingenieurs de reservoir. 

Presentation des figures 

D'autres caracteristiques et avantages de la methode selon 
rinvention, apparaitront a la lecture de la description ci-apres d'un exemple 
non limitatif de realisation, en se referant aux dessins annexes ou : 

- la Fig.1 montre un exemple de schema d'optimisation iterative selon la 
methode; 

- la Fig,2a represente up exemple de champ de permeabilite lognormal 
constituent un modele de reference; 
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- Fig.2b montre les courbes de pression et de sa derivee logarithmique 
issues d'une simulation numerique d'un essai de puits dans le modele de 
reference; 

- la Fig. 3a montre un exemple de champ de permeabilite lognormal associe 
a une premiere realisation du modele; 

- la Fig. 3b montre une comparaison entre les courbes de pression 
associees respectivement a la realisation 1 et au modele de reference; 

- la Fig.4a montre un exemple de champ de permeabilite lognormal associe 
a une deuxieme realisation du modele; 

- la Fig. 4b montre une comparaison entre les courbes de pression 
associees respectivement a la realisation 2 et au modele de reference; 

- les Fig. 5a et 5b representent respectivement les courbes de variation de 
la fonction objectif en fonction du parametre p dans une echelle lineaire 
et dans une echelle logarithmique; 

- la Fig. 6a montre un exemple de champ de permeabilite lognormal associe 
a une realisation contrainte par application de la methode de deformation; 

- la Fig. 6b montre une comparaison entre les courbes de pression de la 
realisation contrainte et du modele de reference. 

DESCRIPTION DETAILLEE DE LA METHODE 

General ites : 

La methode selon ('invention permet de deformer graduellement un 
modele stochastique tout en respectant ses caracteristiques statistiques. 
Cette methode est operationnelle dans un cadre gaussien. Combinee avec 
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une methode ^optimisation, cette methode constitue un outil efficace pour 
calibrer une image du milieu a modeliser a des donnees non lineaires. 

Deformation graduelle aux parametres statistiques fixes : 

On considere un ensemble de realisations possibles R, d'un modele 
stochastique du milieu a modeliser defini par un certain jeu de parametres 
statistiques representatifs (moyenne, variance, variogramme etc). Ces 
differentes realisations R sont obtenues a partir d'un tirage au sort d'un 
germe aleatoire, deux valeurs successives de ce germe pouvant conduire a 
des realisations tres differentes Tune de Tautre. 

La deformation graduelle du modele est obtenue selon la methode par 
une serie de combinaisons iteratives des realisations R^ du modele. Comme 
le montre la Fig.1, la combinaison de base consiste, a partir de deux 
realisations R,. obtenues par tirage. a former une realisation R^, obeissant 
a la relation : 

/?^, =a,./?, (1) 

avec la condition : 

o:r+a2- = l (2) 

Cette realisation combinee se deforme graduellement lorsque Ton 
modifie graduellement les coefficients a, et a. . 

A riteration suivante. on tire une autre realisation R3 et on la combine 
suivant la meme regie avec la realisation combinee precedente R^, 
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Les nouveaux coefficients de la combinaison respectant egalement la 
condition : 

a; +a; = I 

Si le tirage R3. conduit a une realisation combinee qui ne s'accorde 
pas suffisamment avec les mesures, on precede a un ou plusieurs tirages 
successifs jusqu'a I'obtention d'une realisation satisfaisante. 

Le schenna de combinaison precedent peut etre etendu facilement a 
un nombre quelconque de realisations. A la premiere etape de combinaison. 
on tire p realisations differentes et on realise la combinaison suivant la 
relation : 

/? = Xa,/?. (3) 

i 

avec la condition : 

(4). 

A chacune des etapes suivantes, on combine la realisation combinee 
precedente et q autres realisations obtenues par tirage et Ton precede a ieur 
combinaison suivant la meme regie. 

Suivant un mode avantageux de mise en oeuvre, on combine a 
chaque etape deux realisations R,, liees par la relation : 

cos(p7r) + /?, sin(p;r) (5) 
qui respecte la condition (2). 

La minimisation d'une fonction objectif permet de determiner le 
coefficient p . II s'agit d'un probleme d'optimisation mono-parametre. 
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Dans le cas plus general d'utllisation de p+1 realisations, on a un 

problems d'optimlsation des p parametres p, ri = i.2 /;i et les coefficients 

a,. (f= 0.1.2 p) s'ecrivent : 

=sin(p^,7r) 

qui verifient la condition (4). 

Deformation gradueile en modifiant les parametres statistiques: 

Lorsque Ton ne dispose pas de valeurs precises des parametres d'un 
modele stochastique, ces parametres doivent aussi etre identifies par 
rinversion. La methode selon Tinvention permet une deformation gradueile 
d'une realisation d'un modele stochastique en modifiant les parametres du 
variogramme. Pour se faire. nous ecrivons une realisation sous la forme 
suivante : 

ou Y est un bruit blanc (realisation d'un champ gaussien sans correlation) et 
L un operateur determine par le variogramme et permettant de transformer Y 
en une realisation R respectant ce variogramme. Les techniques a mettre en 
oeuvre ici, sont bien connue des specialistes. Une parmi elles est decrite 
par : 

- Oliver, D.S.: "Moving Averages for Gaussian Simulation in two or three 
Dimensions", Math. Geology, Vol.27. No.8, 1995. 
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Nous pouvons ainsi deformer Y de la meme maniere que 
precedemment tout en modifiant Toperateur de variogramme L : 

R = L[Y] ^ L[Y^ cos{ pre) + Y, sin(p7r)] 

ou Y,et Y^sont deux bruits blancs independants I'un de I'autre. 

Deformation graduelle par parties : 

Afin d'ameliorer I'efficacite du calage d'un modele a des donnees non 
lineaires reparties dans le milieu a modeliser, on peut etre amene a diviser le 
milieu en plusieurs parties et effectuer le calage partie par partie. La 
methode selon Tinvention permet d'effectuer des deformations graduelles 
individuelles de differentes parties du milieu a modeliser tout en preservant 
la continuite entre ces parties. Supposons que le milieu soit divise en n 
parties auxquelles sont associees n bruits blancs Y,, Y^, Y„. On peut 
deformer individuellement ces bruits blancs par Talgorithme precedent, et 
puis las rendre correles par Toperateur de correlation L. 



15 



R = l[y] = L 













= L 


y... 


L 



K, cos(p,;r) -h Y^ sin(p,;r) 
K cos(p2;r) -h y; sin(p-,7r) 

y;;cos(p„;r)4-y; sin(p„;r)^ 



ou Y.' et Y." (pour i=1,2,...n) sont deux bruits blancs independants I'un de 
I'autre. 

Extension aux autres modeies: 

La methode de deformation graduelle precedemment decrite peut etre 
20 utilisee dans les cas de modeies transformes a partir d*un modele gaussien 
tels que le modele lognormal et le modele gaussien tronque. 
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Exempi de validation: 

La methode selon ['invention est illustree ici par une etude synthetique 
de la contrainte d'un champ de permeabilite a des donnees de pression 
issues d'un essai de puits. Un channp de permeabilite synthetique (Fig.2a), 
construit selon un modele stochastique lognormal. est utilise comme module 
de reference. Un essai de puits numerique est effectue sur ce champ et it en 
resulte une courbe de pression et sa derivee logarithmlque comme montre la 
figure 2b. II s'agit ici de construire une realisation contrainte par cette courbe 
de pression autre que I'image de reference. 

Les Fig. 3a et 4a montrent deux realisations independantes Tune de 
Tautre du meme modele. On peut constater que les courbes de pression des 
essais de puits sur ces deux realisations sont tres eloignees de celles du 
modele de reference (Fig. 3b et 4b). A partir de ces deux realisations, on a 
construit une chaine de realisations suivant la regie donnee par la formule 
(4). Les Fig. 5a et 5b montrent la variation de la fonction objectif en fonction 
du parametre p (qui varie entre -1 et 1). Les deux realisations initiates 
correspondent respectivement a p = 0 et a p = 0o. On constate que la 
fonction objectif atteint sa valeur minimum pour p = -0.3 environ. En utilisant 
la methode du gradient, il suffit de 3 a 4 iterations pour minimiser cette 
fonction. La Fig. 6 montre la realisation correspondant a p = -0.275 obtenue 
par optimisation et Ton peut constater que le calage a I'essai de puits de 
cette realisation est tres satisfaisant. 
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REVENDICATIONS 

1) Methode pour deformer graduellement une partie au moins de 
realisations d'un modele stochastique de type gaussien ou apparente d*un 

5 milieu heterogene tel qu'une zone souterraine, contraint par un ensemble de 
parametres relatifs a la structure du milieu, caracterisee en ce qu'elle 
comporte le tirage d'un nombre p au moins egal a deux, de realisations 
independantes (Ri) d'une partie au moins du modele choisi du milieu, parmi 
toutes les realisations possibles et au moins une etape de deformation 
10 graduelle du modele en effectuant une combinaison lineaire des p 
realisations initiales independantes, les coefficients (a,) de cette combinaison 
etant tels que la somme de leurs carres est egale a 1 . 

2) Methode selon la revendication 1. caracterisee en ce qu'elle 
comporte plusieurs etapes iteratives de deformation graduelle, avec 

15 combinaison a chaque etape d'une realisation composite obtenue a Tetape 
precedente avec q nouveiles realisations independantes tirees de Tensemble 
des realisations. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte la deformation graduelle d'un modele representatif du milieu tout 

20 en modifiant les parametres statistiques relatifs a la structure du milieu. 

4) Methode selon la revendication 1 ou 2, caracterisee en ce qu'elle 
comporte des deformations graduelles individuelles de differentes parties 
d'un modele representatif du milieu tout en preservant la continuite entre ces 
parties. 
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5) Methode selon Tune des revendications precedentes. caracterisee 
en ce que les coefficients (a) de la combinaison lineaire sont des fonctions 
trigonometriques. 

6) Methode selon Tune des revendications precedentes, caracterisee 
en ce que Ton utilise un modele apparente au type gaussien tel qu*un 
modele lognormal ou gaussien tronque. 
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Fig.2a 
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Fig.2b 
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Fig.6a 
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Fig.6b 
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